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We propose and show that the c-axis transport in high-temperature superconductors is controlled 
by the pseudogap energy and the c-axis resistivity satisfies a universal scaling law in the pseudogap 
phase. We derived approximately a scaling function for the c-axis resistivity and found that it fits 
well with the experimental data of Bi2Sr2CaCu20g+a, Bi2Sr2Ca2Cu30io+* , and YBa2Cu307_a. 
Our works reveals the physical origin of the semiconductor-like behavior of the c-axis resistivity and 
suggests that the c-axis hopping is predominantly coherent. 
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I. INTRODUCTION 

In the pseudogap phase of high-temperature (high-T c ) 
superconductors, the temperature dependence of the c- 
axis resistivity p c is semiconductor-like {dp c /dT < 0), in 
contrast to the in-plane resistivity p ao whose tempera- 
ture dependence is metal-like (dp c /dT > 0). This dra- 
matic difference between p c and p ao is not what one might 
expect within the conventional Fermi liquid theory. It 
has stimulated vast theoretical and experimental investi- 
gations on the interlayer dynamics of high-T c cuprates. 
However, the existing theoretical models based on the no- 
tion of dynamic confinement* 1 ^ or incoherent interlayer 
hoppingS. could not give a natural and unified explana- 
tion for experimental data. 

In high-T c cuprates, the c-axis hopping integral is 
highly anisotropic and depends strongly on the in-plane 
momentai^&l Based on electron structure calculations 4 
and symmetry analysisj&l it was shown that the c-axis 
hopping integral in high-T c cuprates with tetragonal 
symmetry is given by 



t c (k) oc (cosk x — cosky)' 



(1) 



This peculiar in-plane momentum dependence of t c {k) 
results from the hybridization between the bonding O 
2p orbitals and Cu 4s or 3d 3z 2_ r 2 orbitals in each Cu02 
plane. It was confirmed by the angle-resolved photoemis- 
sion spectroscopy (ARPES) measurements^ t c vanishes 
along the two diagonals of the Brillouin zone. These are 
also the directions where the nodes of the d x 2_ y 2 (d)- 
wave superconducting gap or the normal state pseudo- 
gap (it has also the d-wave symmetry2ii£) are located. 
Thus the quasiparticle excitations around the gap nodes 
have almost no contributions to the c-axis hopping. This 
indicates that the c-axis dynamics is governed by the 
quasiparticle excitations around the antinodal points. 



The interplay between t c (k) and the d-wave energy gap 
can affect strongly the c-axis dynamics of electrons. In 
the superconducting state, it leads to, for example, a T 5 
temperature dependent c-axis superfluid density in low 
temperatures^ in contrast to the linear-T behavior of the 
in-plane superfluid density. In a pseudogap normal state, 
since the contribution from the quasiparticle excitations 
around the nodal points is suppressed by the integral (1), 
the interlayer transport would behave as in a gapped sys- 
tem and p c should show a thermally activated behavior. 
On the contrary, the in-plane dynamics is governed by 
low-lying excitations around the gap nodes in low tem- 
peratures. Thus pab would behave similarly as in a gap- 
less system. This explains naturally why the temperature 
dependence of p c is semiconductor-like while p a b is metal- 
like in the pseudogap phase. In making this argument, 
we have implicitly assumed that the c-axis hopping is 
predominantly coherent, 11 We believe this assumption is 
correct, at least in the limit of weak impurity scattering. 
It is supported by the observation of the bilayer splitting 
of the bonding and antibonding bands in the ARPES 
spectra^ the intrinsic T 5 temperature dependence of the 
c-axis superfluid density, as well as the recent angular- 
dependent magnetoresistance measurement^ 

The c-axis dynamics could be also affected by phonon 
or other collective excitations. Along the nodal di- 
rections, a dispersion kink has been observed in the 
spectral function of electrons by the ARPES in both 
superconducting^ and normal states^ This kink is be- 
lieved to result from the electron-phonon coupling. Along 
the antinodal direction, a stronger kink feature was ob- 
served in the single-electron spectrum below T c ^£ This 
stronger kink seems to be correlated with the peak-dip- 
hump lineshape as observed in the ARPES spectra at 
antinodal points^ It is likely to be due to the coupling 
of electrons with a magnetic resonance mode. It could 
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also be explained by the coupling of electrons with a 
Big-phonon mode^£ However, this dispersion kink to- 
gether with the peak-dip-hump structure in the ARPES 
spectra disappears above T c . It suggests that the cou- 
pling between electrons and phonons or magnetic reso- 
nance modes is strong around the gap nodes but very 
weak along the antinodal directions in the normal state. 
Thus we believe that the coupling between electrons and 
phonons or magnetic resonance modes should have very 
weak or negligible effect on c-axis transport properties in 
the normal state. 



II. SCALING HYPOTHESIS 

The above discussion indicates that in the normal state 
the pseudogap is the only energy scale governing the c- 
axis hopping around the antinodal regions. If the c-axis 
hopping is predominately coherent, this would suggest 
that the temperature dependence of p c should satisfy a 
single-parameter scaling law 



T 

p c {T) ota c g[ — 



(2) 



where g(x) is a scaling function, a c is a doping-dependent 
coefficient, and A is the maximal value of the pseudogap. 
The detailed form of g(x) is unknown. However, as will 
be shown below, an approximate but accurate expression 
for g(x) can be obtained by simple theoretical analysis. 

p c is determined by the interlayer hopping integral and 
the scattering mechanism. It can be obtained by evaluat- 
ing the current-current correlation function in the linear 
response theory. In the highly anisotropic cuprate su- 
perconductors, the c-axis hopping can be considered as 
perturbation. In this case, the c-axis conductivity a c is 
given by 



a c (T) oc 



E 

k 



dujv 2 c {k)A 2 {k,uj 



doj 



(3) 



where v c (k) oc t c (k) is the c-axis velocity of electrons and 
f(to) is the Fermi function. A{k,w) is the spectral func- 
tion of pseudogapped electrons. At present, an accurate 
microscopic description of A(k,Lu) is not available since 
the mechanism of pseudogap is still unclear. However, 
phenomenologically one can assume A(k,u>) to have the 
form 



f 



A{k,uu) = - 



TT(0J-E k f+T^ 



(4) 



where Eu = ±y 'e A k + A 2 cos 2 (2</>fc) is the energy disper- 
sion of pseudogapped electrons, <j)k = arctan(fcj, / k x ) and 
Sk is the usual normal state energy measured from the 
Fermi level. T is the linewidth, proportional to the scat- 
tering rate. In the normal state above the pseudogap 
phase, p a b varies linearly with T. It suggests that T scales 
linearly with T. This linear T dependent T is a basic as- 
sumption made in the marginal Fermi liquid thcoryi^. It 



gives a natural account of many experimental results, in- 
cluding the linear in-plane resistivity. Here we will also 
adopt this assumption. Substituting Eq. J3J into Eq. Q), 
it can be shown that to the leading order approximation 
in r, a c is given by 



a c (T) oc -i J duN c {uj) 



dui 



(5) 



where N c (lu) is an effective measure of the c-axis tunnel- 
ing probability of electrons 



(6) 



and V is the system volume. N c (lu) reduces to the density 
of states N{ui) of d-wave gapped electrons if t c (k) in Eq. 
© is set to 1. Fi gure^ compares the energy dependence 
of N c (oj) with N(u>). It shows clearly that the c-axis hop- 
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FIG. 1: Energy dependence of the effective c-axis tunneling 
density of states N c (ui). The density of states of the d-wave 
gapped quasiparticles N(ui) is also shown for comparison. 
N(w) is defined by Eq. © with t c (k) = 1. 

ping of low-energy quasiparticles is strongly suppressed 
by t c {k) and the low energy N c (oj) is vanishingly small 
compared with that at lu ~ A. Thus if the temperature 
considered is not too low compared with A, the contribu- 
tion of low energy electrons to the c-axis hopping can be 
ignored. This has motivated us to replace approximately 
N c (uj) by a step function, 



N c (uj) ~ Nq6(lu — A), 



(7) 



where No is an average value of N c (ui) at lu > A. Here A 
should be slightly smaller than the true maximum pseu- 
dogap since the contribution from the states with lu < A 
should be effectively included in the above approxima- 
tion. Apparently this is a crude approximation, but it 
allows us to get an analytical expression for a c . By sub- 
stituting J7J) into |J5J|, we find that p c is approximately 
given by 



p c (T)=a-\T) 



— r — exp , 
A y \ T 



(8) 
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in the limit T c < T < A. We have tested this for- 
mula with numerical calculations without taking the 
step-function approximation for N c (lo). We find that Eq. 
ijHJ does fit the numerical curves very well in the above 
limit. Though this formula is derived within a limited 
temperature regime, it captures qualitatively the main 
features of p c in the whole temperature range in the nor- 
mal state. In the low temperature limit T -C A, Eq. 
ijSJl is thermally activated. In the high temperature limit 
T » A, p c varies linearly with T. Both agree with the 
experimental observations. 

III. COMPARISON WITH EXPERIMENTS 

To test our single-parameter scaling conjecture, we 
have analyzed the experimental data of p c pub- 



lished for Bi 2 Sr 2 CaCu 2 8+( 5 (Bi2212) by Watanabe et 
al.) 19 i 20 by Chen et al.jSI and by Guira et al.jSS for 
Bi 2 Sr 2 Ca 2 Cu30io + 5 (Bi2223) by Fujii et al.^i and for 
YBa 2 Cu30 7 _«5 (Y123) by Yan et al. 24 and by Babic et 
al.m^ We find that all these experimental data of p c can 
be well fitted by Eq. JSJ nearly in the entire temper- 
ature range. Fig. IHa)-|2fc) show the scaling behavior 
of p c for Bi2212, Bi2223, and Y123, respectively. For 
Bi2212, we only shoed the data published by Watanabe 
et ali 1 ^ for clarity. The experimental data of Bi2212 by 
the same group22i can be also scaled onto the same curve. 
The experimental data of Bi2212 published by other two 
groups 21,22 are showed in Fig. |3| The fitting parameters 
A and a c for Bi2212, Bi2223, and Y123 at different dop- 
ing levels are given in Table I-V. The doping-dependence 
of A and comparison with the ARPES data will be dis- 
cussed later. 




FIG. 2: Scaling behaviors of the c-axis and the in-plane resistivities, (a)-(c), the scaling curves of p c /a c versus T/A for the 
experimental data published by Watanabe et al^- for Bi2212, by Fujii et all— for Bi2223, and by Yan et all— (<5 = 0.3) and by 
Babic et all— — for Y123, respectively. The parameters a c and A are determined by fitting the experimental data with Eq. @ 
well above T c . The dashed lines in (a)-(c) denote the scaling function g(x) = a;exp(l/a;). (d)-(f), the corresponding normalized 
in-plane resistivities [Pab(T)—p a b(0)]/ (a a i>T) as a function of T/A is shown for Bi2212, Bi2223, and Y123, respectively. Different 
curves are shifted vertically from each other for clarity. Both a a b and p a b(0) can be determined by fitting the experimental 
data in the linear regime of p a f, with Eq. I1U1 . In this linear regime, the ratio [p a b(T) — p a b(0)]/(ot a bT) is equal to 1 within 
measurement errors. However, below a doping dependent temperature T*, which is commonly defined as the onset temperature 
of the pseudogap, p a b begins to deviate from this linear-T behavior and [p a b(T) — p a b(0)]/ (a a bT) drops quickly with decreasing 
temperature, as marked roughly by the arrows in (d)-(f). 
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The measurement data deviate slightly from the scal- 
ing curve near T c . This is due to the superconducting 
fluctuations^. This deviation is not unexpected since in 
this case the single-parameter scaling law (2) should be 
modified to include the contributions from the supercon- 
ducting fluctuations. The above analysis shows that the 
scaling behavior of p c is universal and the scaling func- 
tion g(x) is approximately given by 

g(x) = x exp (9) 

independent of the doping concentration as well as 
the chemical structures for these multilayer compounds. 
Thus there is only one energy scale characterizing the 
c-axis dynamics in the pseudogap phase. It suggests 
that the pseudogap is the only energy scale governing 
the quasiparticle excitations at the antinodal points. 



Bi 2 Sr 2 CaCu 2 8+g 




T/A 



TABLE I: The fitting parameters A, a c , a a b, and p a b(0) for 
Bi2212 published in Ref. (the first 7 rows) and Ref. |H 
(the rest 4 rows). The units of a a t and p a b(0) are 10~ 6 Qcm/ K 
and 10 _4 Slcm, respectively. 



s 


T c (K) 


A (K) 


a c (f2cm) 


OLab 


Pab(0) 


0.2135 


71 


416.6 


2.74 


2.60 


1.3 


0.216 


77 


397.4 


2.56 






0.22 


83 


367.3 


2.31 


1.78 


1.1 


0.23 


86 


344.5 


1.95 






0.26 


88 


232.9 


1.07 


1.03 


0.22 


0.27 


85 


207.1 


0.82 


0.78 


0.18 


0.28 


79 


143.3 


0.56 






0.25 


90 


255 


1.14 






0.26 


86 


217 


0.91 






0.27 


83 


184.7 


0.67 






0.28 


79 


150.1 


0.34 






TABLE II: The fittingparameters A and a c 


for Bi2212 pub- 


lished by Chen et al. 
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FIG. 3: The scaled c-axis resistivity for Bi2212 published by 
Chen et al (a)HH and by Sarti et al (b)H2. 



It is interesting to compare the above result with the 
temperature dependence of p a b- In high temperatures, 
p a b shows linear temperature dependence 



Pab{T) = a ab T + Pafe(O), 



(10) 



where a a b is the slope of the linear resistivity and p a b(0) 
is the zero temperature resistivity extrapolated from high 
temperature data of p ab . Fig. Eld)-(f) show the exper- 
imental data of [p a b(T) — p a fc(0)]/(a a bT) as a function 
of T/A for Bi2212, Bi2223, and Y123, respectively. The 
slope a ab and p ab (0) for Bi2212, Bi2223, and Y123 are 
given in Table I, IV, and V, respectively. Here A is the 
pseudogap value determined from p c . An intriguing re- 
sult revealed by Fig. |2Jd)-(f) is that T* (arrows in the 
figure) is approximately proportional to A, i.e., T* cx A. 
This result is reminiscent of the relationship between the 
superconducting transition temperature T c and the su- 
perconducting energy gap Ao. For a BCS mean-field 
d-wave superconductor, T c /Ao ~ 0.47. Here the ratio 
T*/A is about 0.55 for Bi2212, 0.66 for Bi2223, and 0.71 
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TABLE III: The fittin g pa rameters A and a c for Bi2212 pub- 
lished by Giura et al. |22|. 



5 


Tc (K) 


A (K) 


Qc(ficm) 


0.232 


87.4 


321.89 


1.99 


0.240 


88.1 


281.62 


1.64 


0.253 


92.0 


263.84 


1.17 


0.260 


87.0 


194.8 


0.69 


0.265 


85.9 


192.96 


0.54 


0.267 


86.2 


191.69 


0.34 


0.271 


87.6 


177.88 


0.30 



TABLE IV: The fitting parameters A,a c ,a a b, and p a b(0) for 
Bi2223 published in Ref. jSJj. The units of a a b and p ab {0) 
are W~ 6 Qcm/ K and 10~ 5 f2cm, respectively. 



Samples 


T c (K) 


A(K) 


a c (r2cm) 


a a b 


Pab(0) 


a 


87 


403.29 


7.52 






b 


93 


378.74 


6.18 


2.34 


-2.21 


c 


98 


361.74 


5.38 






d 


104 


330.17 


4.15 


2.21 


-4.45 


e 


106 


300.65 


3.4 


2.06 


-5.24 


f 


108 


272.72 


2.78 


1.85 


-5.81 


g 


108 


245.46 


2.28 


1.77 


-6.02 


h 


108 


221.42 


1.85 


1.57 


-4.77 


i 


108 


212.85 


1.03 







for Y123. These values are larger than the BCS value. 
This is not unexpected since, as pointed out above, A 
obtained by Eq. JSJ is smaller than the true pseudogap. 

Figure 0] shows the doping-dependence of A, T*, and 
T c for Bi2212. For comparison, the ARPES data pub- 
lished by Campuzano et ali2£ for A and T* are also shown 
in the figure. The data of T* obtained agree well with 
the ARPES data. Our results of A give a lower bound 
for the maximum pseudogap as expected. They are less 
fluctuating than the ARPES data. It suggests that p c is 
a good probe for the pseudogap. 

The above comparison shows that p c satisfies a simple 
single-parameter scaling law. We believe that the scal- 
ing function g(x) is universal and given approximately by 
Eq. © for all high-T c cuprates whose dominant c-axis 
hopping integral between two neighboring CuC>2 layers 
is coherent and has the in-plane momentum dependence 
given by Eq. Q . This includes all multilayer compounds 
and some single-layer compounds in which Cu atoms 
in the two neighboring unit cells lie collinearly along 
the c-axis, such as HgBa2Cu04+5- However, for other 
single-layer compounds, for example Bi2Si'2CuOg+5 and 
La2- 2; Sr :I ;Cu04, Cu atoms of two adjacent Cu02 planes 
do not lie collinearly along the c-axis. In this case, the 
c-axis hopping integral has the form£iS& 



TABLE V: The fitting parameters A, a c , a a b, and p a b(0) for 
Y123 published in Ref. HH (8 = 0.3) and Ref. El (rest). 

S T c (K) A (K) a c (m£lcm) a ab (p£lcm/ K) p ab (0)(p£lcm) 



0.35 61.4 317.49 10.41 1.36 

0.21 83.7 221.87 3.83 0.86 

0.15 91.0 194.25 2.48 0.84 

0.30 65.0 318.21 105.01 



1.45 
4.77 
4.84 




t c (k) oc cos cos ^ (cos k x 



cos k y Y 



(11) 



0.25 



FIG. 4: Doping dependence of A (Q: Re fs. Il9l20t Q: Ref. 
ED; ® Refjj^), T* (A: Refs. Il9l20ll2 ill22ft . and T c (•: Refs. 
IT9II20II21II22I) for Bi2212. T* is extracted directly from the 
experimental data for the four samples shown in Fig. Hid). 
For all other samples, T* is determined from A by assuming 
T*/A ~ 0.55. The doping concentration p is determined from 
T c using the empirical formula T c /T c>max = 1 — 82.6(p — 0.16) 2 
(dashed line) , T c , max is the maximum superconducting transi- 
tion temperature. For comparison, the ARPES measurement 
data& of A (+) and T* (x) are also shown. The dotted lines 
are guide to eyes. 



It vanishes along both nodal and antinodal directions. 
Therefore the c-axis hopping is dominated by the quasi- 
particle excitations in regions between the nodal and the 
antinodal points. This means that the effective energy 
scale controlling the c-axis dynamics should be the pseu- 
dogap at the momentum on the large Fermi surface con- 
tour where t c (k) takes the maximum, rather than the 
pseudogap itself. In this case, the scaling law Eq. (8) is 
still valid in high temperatures. However, in low temper- 
atures, Eq. (8) needs to be modified since the contribu- 
tion of nodal quasiparticles to the c-axis transport is no 
longer negligible. 



IV. SUMMARY 

The interplay between the anisotropic c-axis hopping 
integral and the pseudogap has profound consequence on 
the c-axis dynamics. It leads to a universal scaling behav- 
ior for p c in the pseudogap phase. The excellent agree- 
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ment between our theoretical analysis and experimental 
data suggests that the pseudogap is the only energy scale 
governing the quasiparticle excitations around the antin- 
odal points in the normal state of high-T c materials. 

Our analysis for the scaling behavior of p c is valid gen- 
erally, independent of the mechanism of normal state 
pseudogap. We believe that it can be generalized and 
applied to other c-axis transport quantities. Before clos- 
ing the paper, we would like to make a general conjecture: 
if F c e is the electronic contribution to a c-axis transport 
coefficient, then F ce should satisfy a single-parameter 
scaling law 

F c , e (T)=/? c / c (£) (12) 

at least at not too low temperatures in the pseudogap 
phase. Here (3 C is a doping-dependent coefficient and 



f c (x) is a corresponding scaling function. A thorough ex- 
amination of this conjecture, no matter being approved 
or disapproved, would shed light on further understand- 
ing of properties of quasiparticle excitations around the 
antinodal points. This will help us to understand more 
about the mechanism of the high-T c superconductiv- 
ity since the superconducting pairing is strongest at the 
antinodal points. 
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